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We study the ground state magnetic properties of ferromagnetic spinor Bose- Einstein condensates 
confined in a deep optical lattices. In the Mott insulator regime, the "mini-condensates" at each 
lattice site behave as mesoscopic spin magnets that can interact with neighboring sites through both 
the static magnetic dipolar interaction and the light-induced dipolar interaction. We show that such 
an array of spin magnets can undergo a ferromagnetic or anti-ferromagnetic phase transition under 
the magnetic dipolar interaction depending on the dimension of the confining optical lattice. The 
ground-state spin configurations and related magnetic properties are investigated in detail. 
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I. INTRODUCTION 

The interaction between quantum degenerate atomic 
gases and optical fields is a cornerstone of modern atomic 
physics and quantum optics. In early experiments on 
Bose-Einstein condensation, light fields were applied pri- 
marily for the capture and precooling of atoms, preced- 
ing the last stage of evaporative cooling in a magnetic 
trap. They were also used to obtain dramatic images of 
the condensates (BEC), and to launch solitons |l| and 
vortices ^ in condensates. This was soon followed by 
applications such as the trapping of condensates in opti- 
cal dipole traps, and the demonstrations of matter- wave 
superradiance ^ Q| and of coherent matter- wave amplifi- 
cation H . More recently, optical dipole traps have been 
employed for the all-optical realization of BEC and of 
quantum-degenerate Fermi gases [|[ . 

Optical lattices, formed by counterpropagating laser 
beams in one, two and three dimensions, were originally 
used in polarization gradient cooling and sub-recoil cool- 
ing experiments at the single-atom level. They rapidly 
found further applications in the manipulation of BECs, 
first in the demonstration of a "mode-locked" atom laser 
and the observation of Joscphson tunneling between lat- 
tice wells H, and subsequently in the transport and ac- 
celeration of condensates |J. More recently, they have 
lead to the demonstration of the superfluid-Mott insu- 
lator transition |ll), and of the collapse and revival 
of the condensate wave function ]l2| . In the near fu- 
ture, they may also prove useful in the realization of 
bright atomic solitons relying on negative effective atomic 
masses in periodic potentials 

In contrast to magnetic traps, which only capture 
atoms in weak-field seeking states, optical traps function 
for all hyperfine sublevels of the alkali electronic ground 
states. This presents considerable advantages, in partic- 
ular in the study of spinor condensates such as sodium 
and rubidium. The first study of the magnetic proper- 
ties of spinor condensates were carried out by Ketterle 
and coworkers, who investigated the existence of coex- 
isting spin domains in 23 Na, an "anti-ferromagnetic", or 
"polar" condensate [0. 

Recent experimental and theoretical studies have es- 



tablished that in contrast to 23 Na, 87 Rb is expected to be 
ferromagnetic at zero temperature. That is, the expecta- 
tion value of its total spin F is finite, (F) ^ [||[ |J, 0. 
As a result, an ensemble of condensates placed at the po- 
tential minima of an optical lattice would act as meso- 
scopic magnets, much like large spins on a crystalline 
lattice. In the absence of external fields and long range 
site-to-site interactions, these magnets would have ran- 
dom orientations. 

The situation is changed in the presence of interac- 
tions between neighboring lattice sites. It is known that 
in the case of spins on a crystal lattice, the dominant 
source of coupling is the quantum-mechanical exchange 
interaction. We recall that 19th century physics failed in 
its attempts to explain ferromagnetism in terms of the 
magnetic dipole-dipole interaction, and it is Heisenberg 
who first introduced the exchange force to explain this 
effect 
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In the present case, though, the overlap 
between neighboring condensate wave functions is neg- 
ligible for deep enough lattice wells — the Mott insula- 
tor state — and so is the exchange interaction. Instead, 
the individual mesoscopic magnets are coupled by the 
magnetic (and possibly also the optical) dipole-dipole in- 
teraction. Because of the large number N of atoms at 
each lattice site, this interaction is no longer negligible, 
despite the large distance, of the order of half an optical 
wavelength, between sites. As such, the present situation 
is in some sense a return to 19th century physics. The 
goal of this paper is to discuss several aspects of the spin 
and magnetic properties of such lattice systems in one- 
and two-dimensions. 

The remainder of the paper is organized as follows. 
Section II briefly reviews the theory of spinor conden- 
sates in general, with special emphasis on the ferromag- 
netic and polar ground states resulting from local spin- 
changing collisions. We then introduce the nonlocal long- 
range magnetic dipole-dipole interaction between con- 
densates at different sites in the optical lattice. Sec- 
tion III briefly reviews previously published results on 
one-dimensional lattices, and discusses the ferromagnetic 
ground state of the full lattice. On this basis, an ex- 
tension of the one-dimensional case to two-dimensional 
lattices are analyzed in section IV. The ground state of 
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the system is determined numerically using a genetic 
algorithm that is discussed in some detail. We show 
that in that case, the ground state is normally anti- 
ferromagnetic. Edge effects are also briefly addressed. 
Finally, Section V is a conclusion and outlook. 



II. MODEL 

The dynamics of spinor condensates trapped in optical 
lattices is primarily governed by three types of two-body 
interactions: spin-changing collisions, magnetic dipole- 
dipole interactions, and light-induced dipole-dipole inter- 
actions. For an optical lattice created by blue-detuned 
laser beams, the atoms are trapped in the dark-field 
nodes of the lattice and the light-induced dipole-dipole 
interaction can be neglected Q. In this paper, we focus 
on this case. As a preparation for sections III, IV and 
V, we first discuss the interatomic interactions in some 
detail. 



A. Spin-changing collisions 

In second-quantized form, the Hamiltonian describing 
a system of spin / = 1 bosons subject to local spin- 
changing collisions is p5J ^3|, Q 



[ d 3 r^(r)Vj(r)^(r)^(r) 



(1) 



+ f E /rf 3 ^(r)^(r)F a , M -F^,^(r)^(r), 
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where tp a (r) is the field annihilation operator for an atom 
in the hyperfine state |/ = l,m/ = a), a = —1, 0, 1, U(r) 
is a potential produced by an optical dipole trap and as- 
sumed to be the same for all hyperfine states, and M is 
the mass of the atoms. F is the vector operator for the 
hyperfine spin of an atom, with components represented 
by 3 x 3 matrices in the \ f — {,mf — a) subspace. For 
ultracold bosons, only s-wave collisions with total hyper- 
fine spin of F — 0, 2 are allowed, and 
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where ao and a 2 are the s-wave scattering lengths for 
collisions in the F = and F = 2 channel, respectively. 

The ground state properties of spinor condensates sub- 
ject to these local spin-changing collisions have been de- 
termined by introducing the components a (r) of the 



spinor condensate wave function in the mean-field ap- 
proximation, 



</>a(r) = (0* (r)) = y/n{r)C a (r), 



(2) 



where n(r) is the local atomic density and Ca( r ) a nor ~ 
malized spinor, and minimizing the energy functional 
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In this expression, /i is the chemical potential and the 
averaged single-atom spin angular momentum is 
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For c 2 > 0, the energy E is minimized by (F(r)} = 0, 
and the spinor condensate is in an "anti-ferromagnetic" , 
or "polar" state. This is the case for 23 Na condensates, in 
which case a 2 — ao — 5 a.u. Ketterle and coworkers have 
studied this situation in great detail [l4| . In particular 
they have obtained spin-domain diagrams and studied 
experimentally the miscibility of these domains in the 
presence of external fields. 

For c 2 < 0, in contrast, the energy E is minimized 
by making (F(r)) 2 = 1. As discussed in Ref. 22 , the 
direction of the spin is 



(F(r)} = cos (3qz + sin (3q x (cosaox + sinaoy), 



(4) 



where ao and f3 are Euler angles. All possible ori- 
entations {a ,P ) are possible and lead to the same 
ground-state energy E. Recent theoretical calculations 
by Klausen et al. predict that for spin-1 87 Rb, the 
scattering lengths ao and a 2 are almost equal, but with 
ao > a 2 , with a difference of the order of 0.3 to 2.7 a.u. 
& 

Consider, then, an 87 Rb condensate trapped on an op- 
tical lattice with wells deep enough that its ground state 
is the Mott-insulator state, i.e., there is no global phase 
of the condensate over many lattice sites [[□] . Each lat- 
tice site is therefore the location of a "mini-condensate," 
which can contain as many as several thousands atoms 
in one-dimensional lattices, and several hundreds in 2- 
D lattices. In the absence of external fields and long- 
range site-to-site interactions, these condensates can be 
thought of as independent magnets, whose spin vectors 
point in random directions, with no spin correlations be- 
tween sites. This situation is similar to the spin lattices 
familiar from the study of magnetism, with two differ- 
ences. First, the quantum mechanical exchange interac- 
tion, which is at the core of magnetism, is completely neg- 
ligible in the present situation. This is because neighbor- 
ing sites on an optical lattice are at least one half optical 
wavelength apart. For deep lattice wells, the center-of- 
mass wave functions for the individual mini-condensate 
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— essentially the ground state Wannier wave functions at 
the individual sites — do not have any significant overlap. 
Second, the magnetic dipolar coupling, which is normally 
negligible and leads to the prediction of Curie tempera- 
tures several orders of magnitude lower than actually ob- 
served in solid state magnetic materials, is now the dom- 
inant interaction, due to the large number TV of atoms 
at each lattice site. This leads to an N 2 enhancement 
factor, as we now show. 



B. Magnetic dipole-dipole interaction 

In order to describe the magnetic dipolar interaction 
between mini-condensates at lattice sites i and j, we as- 
sume that the condensates at each site can be treated 
independently, and have the same spatial form, which is 
also independent of the spin state of the atoms. Specifi- 
cally, we decompose the Schrodinger field operator as 



V>(r) 



with 



E 
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In this expansion, which goes beyond the mean-field ap- 
proximation of Eq. (||), rj is the coordinate of the i-th 
lattice site, a a (i) and aj,,(i) are bosonic annihilation and 
creation operators for atoms in the hyperfinc state a at 
site i, and </>j(r) = tf>(r — r ; ) is the ground state wave 
function of the mini-condensate at that site, normalized 
to unity. For ag ~ a2, it is approximately given by the 
solution of the stationary Gross-Pitaevskii equation 
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where 
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is the total number of atoms at site i and we assume that 
all sites have the same number of atoms. 

The magnetic dipole-dipole interaction between the 
mini-condensates at sites i and j is given by |25|] 



= £ J d 3 r j dV |0(r - r ; )|W - rj)| 
Mi-Mj 3(/Xi-(r-r'))(Mj-(r-r')) 



r r\3 



./|5 



where fio is the vacuum permeability and fli is the mag- 
netic dipole moment at site i. In second-quantized form, 
it is given explicitly by 

Mi = ■y B ^2al t (i)F a . l3 ap{i) = 7 B S ; , 



where 73 = qf^b is the gyromagnetic ratio and we rec- 
ognize that S; is the angular momentum operator for the 
condensate at site i. We remark that for a given site, the 
expectation value of /2i is 

(Mi) = lB^2(ai(i)F a ^a/3(i)) 

a,0 

where (Pi) is the single-atom magnetization at the site, 
see Eq. (g). 

Summarizing, then, the Hamiltonian describing the 
spinor "mini-condensates" in the optical lattice, subject 
to spin-changing collisions and to an inter-site magnetic 
dipolar interaction has the spin-spin coupling form 
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and the tensor Ay is defined by 
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We have also introduced an external magnetic field B cxt 
for future use. In the limit of tight confinement, the 
condensate wave functions at each lattice site can be ap- 
proximated by 

|^(r-ri)| 2 »5(r-n). 
In this limit we have 

, _ IB Mo 
j ~ 47r|r y |3' 

and the tensor Ay becomes 

^ij = ^ij*fj! 

where ry = r ; — rj and fy = ry/|ry|. 



III. FERROMAGNETISM IN A ID OPTICAL 
LATTICE 

In this section, we study the magnetic properties and 
spin dynamics of spinor condensates in a ID optical lat- 
tice. More specifically, we consider a blue-detuned opti- 
cal lattice where the mini-condensates are trapped at the 
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standing wave nodes. In this case, the light-induced dipo- 
lar interaction can be ignored and the mini-condensates 
only interact via the magnetic dipolar interaction. With- 
out loss of generality, we assume that the axis of the 
lattice is along the z direction, which we also choose as 
the quantization axis. Hence the total Hamiltonian (|^) 
reduces to 
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We assume that the magnetic field B ext is of the form 



B, 



= B Z 



B x x, 



where B z z is an applied field and B x i is an effective 
magnetic field that accounts for all possible effects from 
the experimental environment. While this field can have 
any possible orientation, we take it to be transverse and 
along x without loss of generality, since any longitudinal 
component can be included in B z . 

Furthermore, we consider an infinitely long lattice so 
that boundary effects can be ignored. The hamiltonian 
describing the spin S of a generic site i reads then 



h = \' a S 2 - 1B S 



B~ 



(8) 



We now proceed to determine the ground state of the 
single-site Hamiltonian (|J) in the mean-field — or Weiss 
molecular field — approximation 19 . It consists in re- 
placing the operators S", a = x,y,z, by their ground- 
state expectation value 



(St) 



Nm a , 



(9) 



which is assumed to be the same for all sites. We re- 
mark that m z is nothing but the difference in population 
of the Zeeman sublevels of magnetic quantum numbers 
±1. Replacing by Nm a allows us to approximate the 
Hamiltonian (||) by 

/ lmf -A^S 2 - 7B S-B cff , (10) 

where we have introduced the effective magnetic field 

B c ff = (B z + 2Am z )i + (B x - Am x )x - Am y y, 

and 

A = N"£\ ij . 



In the case of 87 Rb, the individual spinor condensates 
at the lattice sites are ferromagnetic, X' a < 0. In that 
case, the ground state of the mean-field Hamiltonian ( |l0| ) 
must correspond to a situation where the condensate at 
the site i under consideration must be aligned along B e a 
and takes its maximum possible value N. That is, the 
ground state of the mean-field Hamiltonian (|l0|) is simply 



\GS) = \N,N)i 



(11) 



where the first number denotes the total angular momen- 
tum and the second its component along the direction of 
B c ff. Note that \GS) represents a spin coherent state 
in the basis of \S,S Z ). The fact that the ground state 
magnetic dipole moment of each lattice site is N times 
that of an individual atom results in a significant mag- 
netic dipole-dipole interaction even for lattice points sep- 
arated by hundreds of nanometers. This feature, which 
can be interpreted as a signature of Bose enhancement, 
is in stark contrast with usual ferromagnetism, where the 
magnetic interaction is negligible compared to exchange 
and where the use of fermions is essential. 

The mean-field ground state of Eq. ([ll]) allows us to 
calculate the magnetization m x ,y,z- One finds readily 

m a = -^{GS\S?\GS)=cos6 a , 

where 9 a is the angle between B e ff and the a-axis. In the 
absence of externally applied field, B z = 0, this gives 
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where B = ^(2Am z ) 2 + (B x - Am x ) 2 + (Am y ) 2 nor- 
malizes the magnetization vector to unity . 

Since B > 0, the third of these equations implies 
that Ttiy = 0. With m 2 + m 2 . = 1 and the condition 
2 A = B, which follows directly from the equation for m z , 
we find further that for B x > 3A, the unique solution is 
m z = m y = 0, m x = 1. That is, the lattice of conden- 
sates is magnetically polarized along the environmental 
magnetic field B x . For B x < 3A, in contrast, there are 
two coexisting sets of solutions: i) m z = m y — and 
m x = 1; and ii) m z — ±yl — (-Ba;/3A) 2 , m y = and 
rrix — B X /3A. It is easily seen that the state associated 
with the latter solutions has the lower energy . Hence 
it corresponds to the true ground state, while solution 1 
represents an unstable equilibrium. 

We have, then, the following situation: As the effective 
magnetic field strength B x is reduced below the critical 
value 3A, the lattice ceases to be polarized along the 
direction of that field. A phase transition occurs, and 
a spontaneous magnetization along the z-direction ap- 
pears, characterized by a finite m z . This phenomenon 
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is reminiscent of conventional ferromagnetism. Indeed, 
our model is analogous to the Ising model with the 
environmental transverse magnetic field B x playing the 
role of temperature. For B x — — corresponding to zero 
temperature in Ising model — the spins at each lattice 
site Si align themselves along the lattice direction, even 
in the absence of longitudinal field. This spontaneous 
spin magnetization diminishes as B x increases, and com- 
pletely vanishes if B x exceeds the critical value 3A — the 
analog of the Curie temperature in the Ising model. We 
note however that the situation at hand exhibits impor- 
tant qualitative differences with the Ising model. For ex- 
ample, no spontaneous magnetization occurs in ID Ising 
model, for any finite temperature. 

We note however that the appearance of a spontaneous 
magnetization does not rely on this condition being ful- 
filled. This point was discussed in Ref.^Tj, which nu- 
merically solved the Hamiltonian (0) without invoking 
the mean-field approximation for a two-well system and 
showed how the situation rapidly approaches the mean- 
field results as N increases. 



IV. ANTI-FERROMAGNETIC GROUND STATE 
OF THE 2D LATTICE 

Now we turn our attention to 2D lattices, formed as 
before by blue-detuned lasers. We show that depending 
on the relative magnitude of the lattice constants along 
its two axes, this system exhibits a variety of possible 
ground states, including an anti-ferromagnetic configu- 
ration. 

We consider a rectangular lattice in the (y, z)-plane, 
with primitive lattice vectors a = az and b = by, of 
lengths a and b, in these two directions. We assume as 
before that the number of atoms at each lattice site is 
the same and that the atoms are tightly confined so that 
we can approximate their probability density by a delta 
function at each lattice site, 

|^(r)| 2 = 5(r-r^). 

Here, = i&+jh is the position of the center of the 
lattice site. Under these conditions, the Hamiltonian (||) 
with B oxt = becomes 



H 



where 
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and Yij t ki = r ij ~ r ki = na + fib, with n = i — k and 
m = j — I. 



A. infinite size lattices 

As in the preceding section, we determine the ground 
state of the lattice in the semiclassical limit, ignoring 
spin-spin correlations and replacing the operators Sy- 
with their expectation value with respect to a spin co- 
herent state, 



The semiclassical ground state corresponds to the orien- 
tation of the spin vectors that minimizes the semiclassical 
energy corresponding to the Hamiltonian JT3]). In con- 
trast to the one-dimensional case, it is not obvious from 
inspection of Eq. (|l3|) that all expectation values My 
should be equal. Hence, the determination of the ground 
state for an N x M lattice requires the minimalization 
of the energy with respect to 2NM variables. However, 
in the limit of an infinite lattice the ground state should 
be translationally invariant with respect to displacements 
of the spins by a finite number of lattice constants along 
either axis. We can therefore generalize the mean-field 
ansatz used in the one-dimensional case by assuming that 
the 2D lattice can be decomposed into a finite set {£} of 
interpenetrating periodic sublattices for which all spin 
vectors have the same orientation. 

The positions of the sites of the sublattice I of primitive 
lattice vectors en and are 

rt t ij = ia e + & efi + jb e + b eM 

where a^o and b^.o denote the origin of that lattice. Since 
the interaction between dipole moments that are perpen- 
dicular to the plane of the lattice is repulsive while the 
interaction between dipole moments in the plane of the 
lattice is predominantly attractive, the ground state must 
correspond to spin vectors in the plane {y, z) of the lat- 
tice. Hence the spin vector associated with the sublattice 
I can be written as 

M e = N (cos 9 e y + sin 6 e z) . 

One can gain an intuitive feel for the ground state of 
the system by considering what happens when one lets a 
ID lattice approach an already existing one from infin- 
ity. For concreteness, we take the axes of both lattices to 
be along z. We know from the previous section that for 
large lattice separations, the spins in each lattice will be 
oriented in either the +Z or — z direction with equal prob- 
ability. In effect, each lattice acts like a long bar magnet. 
As the lattices approach each other, though, they start 
to interact via their magnetic dipole moments. The min- 
imization of energy then proceeds in a familiar way: Just 
as two bar magnets placed side by side orient themselves 
so that opposite poles are next to each other, the spins 
of the two ID lattices will arrange their orientation so 
that the spins in one lattice point along +z while in the 
other lattice the spins point along — z. This will remain 
true as long as the lattice separation is much larger than 
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the primitive lattice vector b of the ID lattices, due to 
the 1/r 3 dependence of the magnetic dipole interaction. 
Indeed, in this case the easy axis is the y-axis. 

This argument can easily be generalized to many rows. 
It follows that for a ^> b, rows of spins parallel to the pr- 
axis will alternatively align themselves along the +z and 
— z direction. Similarly, for o « d, the z-axis becomes 
the easy axis and rows of spins parallel to y align them- 
selves alternatively along the +y or — y direction. In 
both cases, though, the ground state is expected to be 
anti- ferromagnetic. 

Even though the magnetic dipole interaction is long 
ranged, it is easy to see that neighboring spins within 
each row interact more strongly than do neighboring 
spins in adjacent rows provided that a = b + e with e 
positive. One therefore expects that the ground state 
will remain anti-ferromagnetic unless e — > 0. In this 
limit, there are clearly two degenerate anti-ferromagnetic 
ground states that are topologically distinct, i.e. that can 
not be related by a simple rotation. Any weighted com- 
bination of these two configurations has the same energy 
and is therefore a new degenerate ground state. Assign- 
ing the weight cos 2 9 to the ground state with all spins 
pointing in the ±y direction and sin 2 9 to the ground 
state with all spins in the ±z, then we find that this sit- 
uation is equivalent to a ground state consisting of four 
interpenetrating sublattices [I = 1, 2, 3, 4), with spin ori- 
entations, 

Mi = AT(cos0y + sin0z), (15a) 
M 2 = N(-cos9y + sin#z), (15b) 
M 3 = N(-cos9y - sinflz), (15c) 
M 4 = N(cos9y - sinflz), (15d) 

and sublattice sites located at 

ri t ij = 2ia+2jb, (16a) 

r 2 , ij = 2za+(2j + l)b, (16b) 

(2i + l)a+(2j + l)b, (16c) 

r 4l y = (2i + l)a+2jb, (16d) 

where i,j — 0, ±1,±2, .... The corresponding lattice 
structure is illustrated in Fig. |l|. For a = b all values 
of 9 are degenerate while for a < b and a > b the ground 
state corresponds to 9 equal to and 9 = ir/2, respec- 
tively. 

The next section discusses the use of a genetic al- 
gorithm to numerically determine of the lattice ground 
state for a finite lattice size. 



B. Finite size lattices 

1. Genetic algorithm 

Genetic algorithms have become a widely used tool 
for solving optimization problems that depend on a large 




0.5 

Z 

FIG. 1: Orientations of the spins on the four interpenetrating 
sublattices for a = b. The lengths are in units of A. 

number of variables p8[ . The basic idea behind genetic 
algorithms is Darwinian natural selection. These algo- 
rithms proceed from an initial set of trial solutions to 
the optimization problem, which can be thought of as in- 
dividuals in a population. The individuals breed, follow- 
ing some prescribed mating rules, to produce offspring, 
which constitute the next generation of individuals. In 
addition, random mutations are also introduced. The 
offspring that produce better solutions to the problem 
survive and are allowed to further breed, while those 
that produce poor solutions are eliminated. Ideally, after 
many generations the algorithm converges to the optimal 
solution(s) to the problem at hand. 

In the specific system at hand, the algorithm starts 
from a large population M of initial lattices, typically 
Af = 512. Most of them have completely random spin 
orientations, but some may have ordered configurations 
based on the ground state of the infinite lattice. At each 
generation, the genetic algorithm performs a combination 
of mutations and breeding steps on the members of the 
population, which we refer to as mutating and mating. 

The mutations modify each member of the population 
to form a second population of N lattices. They can be 
either global and local. Local mutations involve giving 
random rotations to a random percentage of the spins in 
the individual lattices. These rotations are by angles (p 
and 9 about the y and x axes, respectively, where ip and 9 
are normally distributed random numbers with standard 
deviations typically chosen to be n/8. In contrast, the 
global mutations rotate all spin in the lattice by related 
amounts: They either apply the same random rotation 
to all lattice sites, or rotate the spin at each lattice site 
by a slightly different amount determined by its value 
(this is used when investigating the case of equal lattice 
constants, a — b, discussed below). In general, a given in- 
dividual is subjected to both local and global mutations. 

After the mutations are performed, the 2AF individu- 
als are allowed to mate. The mating process randomly 
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picks two individuals using a normally distributed prob- 
ability distribution centered around individuals with the 
lowest energy. This insures that, on average, only those 
individuals with the lowest energies produce offspring. 
Each pair of parents produces two offsprings using one 
of four randomly selected mating techniques: site swap- 
ping, sub-lattice swapping, row and column swapping, 
and row and column rearranging. Site swapping consists 
of swapping a random number of randomly chosen sites 
from the parents. Similarly, sub-lattice swapping consists 
of swapping a randomly sized and positioned sub-lattice 
between the parents. Row and column swapping works 
by randomly picking rows from both parents and form- 
ing one child, and doing the same with columns to form 
a second child. Row and column rearranging uses only a 
single parent to produce a child by randomly rearranging 
its rows or columns. The mating process is repeated Af 
times at each generation to produce a total population 
of 4JV lattices. Of those, only the TV individuals with the 
lowest energy are selected as parents for the next gener- 
ation. 

The genetic algorithm is run until the relative energies 
of the individuals in generation M and M. — 100 differ 
by less than 10~ 7 . 



2. Numerical results 

The ground state of the system determined by the ge- 
netic algorithm is characterized by all spins lying in the 
plane of the lattice, in agreement with the discussion of 
section IV. If the lengths of the primitive lattice vectors a 
and b differ significantly, say, by 10 percent or more, the 
ground state is anti- ferromagnetic. With the exception 
of sites near the lattice boundary, the anti-ferromagnetic 
structure is identical to that predicted based on an infi- 
nite lattice. 

As is to be expected, boundary effects become more 
important, the smaller the lattice. In that case, the 
ground state is characterized by spins orientations near 
the boundaries that deviate from the ±y or ±z direc- 
tions. When a and b are significantly different, these 
boundary effects are manifest only near the corners of 
the lattice, and they lower the ground state energy by 
a very small amount. For example, for a = 0.6A and 
b = 0.5 A, where A is the wavelength of the laser forming 
the lattice in the z-direction, the boundary effects reduce 
the ground-state energy of an 11 x 11 lattice by only 0.1 
percent compared to its infinite lattice value. For larger 
lattices, the boundary effects become even smaller. 

When a = 6, finite size effects are more important in 
determining the spin structure of the ground state. We 
recall that in that case, an infinite lattice possesses an 
infinite number of degenerate ground states characterized 
by the angle 0. Boundary effects break this degeneracy 
and lead to the appearance of a preferred pattern. Fig. ^ 
illustrates the transition from the boundary dominated 
pattern of the a = b situation to the anti-ferromagnetic 
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FIG. 2: The ground state configuration shows the transition 
from the boundary dominated pattern for a — b to the anti- 
ferromagnetic configuration for a ^ b. From top to bottom, 
b = 0.5, 0.505 and 0.6, respectively and a = 0.5 for all figures. 
The lengths are in units of A. 



configuration of a ^ b. As illustrated in Fig. ga, for 
the case of a = b, near boundaries the spins are aligned 
parallel to them. That this should be the case is plausible 
since when going from a situation where a < b to a > b, 
the spin orientation must go from being parallel to the y- 
axis to being parallel to the z-axis. To accommodate the 
orthogonal directions along two adjacent boundaries, the 
angle near the corners changes in such a way that the 
spins at the corner sites make an angle of ir/4 relative 
to the y- and z-axis. This lifts the degeneracy present 
in the infinite lattice. As a result, the spins near the 
center of the finite lattice always take on an orientation 
corresponding to Eqs. ( |i"5| ) and ( |l6| ) with = 7r/4. This 
result holds for all finite-size lattices. Finally, we note 
that the ground-states of finite-size lattices are two-fold 
degenerate, the second ground state being obtained by 
reflections about the y and z axes. 

Figure ^ shows how the spins orients themselves as b 
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FIG. 3: Plot of the deviation angle 6 relative to the y axis for 
the spin at the center of the lattice (solid curve) and a spin 
on left boundary of the lattice (dashed curve) as functions of 
b for fixed a — 0.5. Lengths are in units of A. 

changes for fixed a. As b deviates from a, the spins near 
the center of the lattice quickly become parallel to the 
easy-axis, while the spins near the boundaries become so 
much more slowly. 

V. SUMMARY AND OUTLOOK 

In summary, we have studied the spin configurations 
and magnetic properties of spinor Bose-Einstein conden- 
sates in an optical lattice. In the tight-binding limit, 
the ground state is the Mott-insulator state and the con- 
densed atoms at each lattice site collectively behave as 
a spin magnet. Due to Bose enhancement, the dipole- 
dipole interactions between these spin magnets become 
important and may give rise to a rich variety of phenom- 
ena. We have shown here that the array of spin magnets 
can undergo a ferromagnetic (in the ID case) or anti- 
ferromagnetic (in the 2D case) phase transition under 
the dipolar interaction when external magnetic fields are 
sufficiently weak. Using the same mechanism, it will also 
be possible to create ferrimagnetic lattice systems if one 
can interleave two sets of optical dipole potentials, each 
trapping one species of atoms (or one hyperfine state of 
the same atom) different from the other. 

In the case of a far red-detuned lattice such that the 
spacing between adjacent lattice site exceeds the atomic 
resonant wavelength, the detection of the ground state 
spin structure amounts to detecting populations in the 
individual Zeeman sublevels at each site. This can be 
achieved using a Raman scattering scheme. For example, 
one can shine two light beams, one 7r-polarized and the 
other circularly polarized, onto the system. The absorp- 
tion or gain of the probing light after passing the sample 



is then a measure of the relative population of the hy- 
perfine levels, since it depends upon which of them are 
initially populated. This scheme wouldn't work for a blue 
detuned lattice, though, since in that case the spacing be- 
tween neighboring sites is sub-wavelength. However, the 
long range periodic spin structure, in particular the fer- 
romagnetic and anti-ferroma gne tic ordering, can still be 
detected by Bragg scattering |2!|. Let us take 87 Rb as an 
example. Its ground state is the 5Si/2 state with F = 1. 
For cr + -polarized Bragg probe light (we choose the quan- 
tization axis to be parallel or antiparallel to atomic spins) 
with a frequency close to the F = 1 — > F' = 2 D2 reso- 
nance line , then the ratio of the transition strength (or 
scattering cross section) for atoms in m = —1 and m = 1 
Zeeman sublevel is 1/6. As a result, the Bragg signal de- 
pends on whether one has a ferromagnetic lattice (where 
all the atoms are in either m = — 1 or m = 1 Zeeman sub- 
level) or anti-ferromagnetic lattice (where half the atoms 
are mm — — 1 and the other half are in m = 1 sublevel) . 

In addition to their ground state structure, spinor con- 
densates in an optical lattices also possesses consider- 
able potential for studying other phenomena such as spin 
waves ||^] , macroscopic magnetization tunneling |3l|| , do- 
main wall formation, etc. Future studies will also include 
the dynamical properties of the system. Due to the long- 
range as well as the nonlinear nature of the dipolar inter- 
action, the dynamics of the system should be very rich. 
For instance, given a ground state 2D lattice with primi- 
tive lattice constants a < b where all the spins are aligned 
along the ±y direction, one can suddenly modify the lat- 
tice light so that a > b. Whether and how the spins 
adjust themselves to the new ground state will be an in- 
teresting problem, closely related to the phenomenon of 
spin tunneling |3l| . In addition, these systems may also 
find applications in the field of quantum information and 
computation. We conclude by noting that in addition to 
the Mott insulator limit studied in this paper, the ge- 
netic algorithm that we have developed here might be 
modified to investigate the other limit where tunneling 
between lattice sites becomes significant and the system 
becomes a superfluid J|2| . 
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